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ABSTRACT 

Building on our previous spectroscopic and photometric analysis of the isolated Local Group dwarf 
irregular (dlrr) galaxy WLM, we present a comparison of the metallicities of its RGB stars with respect 
to the well studied Local Group dwarf spheroidal galaxies (dSphs) and Magellanic Clouds. We calcu- 
late a mean metallicity of [Fe/H]= —1.28 ± 0.02, and intrinsic spread in metallicity of a = 0.38 ± 0.04 
dex, similar to the mean and spread observed in the massive dSph Fornax and the Small Magel- 
lanic Cloud. Thus, despite its isolated environment the global metallicity still follows expectations 
for WLM's mass and its global chemical evolution is similar to other nearby luminous dwarf galaxies 
(gas-rich or gas-poor). The data also show a radial gradient in [Fc/H] of d[Fc/H]/dr c = —0.04 ± 0.04 



dex 



which is flatter than that seen in the unbiased and spatially extended surveys of dSphs. 



Comparison of the spatial distribution of [Fe/H] in WLM, the Magellanic Clouds, and a sample of 
Local Group dSphs, shows an apparent dichotomy in the sense that the dlrrs have statistically flatter 
radial [Fe/H] gradients than the low angular momentum dSphs. The correlation between angular 
momentum and radial metallicity gradient is further supported when considering the Local Group 
dEs. This chemodynamic relationship offers a new and useful constraint for environment driven dwarf 
galaxy evolution models in the Local Group. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: dwarf — galaxies: Local 
Group — galaxies: individual (WLM) 
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1. INTRODUCTION 

Dwarf galaxies in the Local Group offer a strong test 
for mass assembly theories in ACDM cosmologies, as hi- 
erarchical merging of protogalactic fragments of similar 
stellar mass (10 5 — 10 9 M ) are expected to be one chan- 
nel for formation of larger disk galaxies like our own 
Milky Way (MW) or Andromeda (M31) (e.g., |Navarro 
et al.|1997[|Moore et al|1999||Madau et al.|2001[ ). Char- 
actenzmg the physical properties of the building blocks 
of larger galaxies requires an understanding of the global 
properties of dwarf galaxies such as mass, size, angu- 
lar momentum, chemistry and luminosity, to observe if 
the merging fragments are consistent with the properties 
of the disk and halo of the large spirals (however it is 
not certain how sim ilar the surviving dwarfs are to those 
merging fragments; Font et aL]|2006[ ). 

Studies of dwarf galaxies can also shed light on how 
baryons populate dark matter halos at the faint end of 
the galaxy luminosity function, and in turn offer a chance 
to study the star format ion and galaxy assembly at low 
metallicities and masses (iBullock & Johnston 20051. If 
low star formation efficie ncy m d warf galaxies is due to 
internal feedback effects dBrooks et al.|| 2007t, or H2 reg- 

DekV 



also provide insight towards what produces the range 
of morphologies (dwarf spheroidals, transition dwarfs, 
dwarf irregulars) seen in the Local Group. Quantitative 
comparisons of the chemical and kinematic signatures of 
dwarf galaxies can also constrain evolutio nary connec- 
tions (e.g., tidal transformation scenarios; |Mayer et al.| 
2006| ) between the two morphological classes of dwarf 



ulated st^ovma,tion(l^^^^Wmiem^ wfll galaxies the recent study o flWeisz et al.|fl2011a| showed 



galaxies. 

Detailed signatures of environmental or internal evolu- 
tionary mechanisms can be traced by the evolved stel- 
lar populations (RGB stars) in dwarf galaxies. With 
spectroscopic observations of these long lived (> 1 Gyr) 
stars, questions on whether star formation and chemi- 
cal enrichment proceed with the same efficiency in iso- 
lated and tidally perturbed galaxies can be addressed 
through analysis of the stellar age-metallicity relation 
and star formation history (SFH). While the present day 
appearance and morphology-density relations exhibited 
by dwarf galaxies i n the Local Group and other groups 
( Weisz et al.|2011b ) indicate that they are susceptible to 
environmental processing, the details of how environment 
has influenced them over a Hubble time are still difficult 
to ascertain. For example, in a sample of distan t dwarf 



the age and metallicity of a galaxy's stellar population 
change substantially? Understanding this question can 

rleaman@iac.es 



that the SFH of dlrrs and dSphs are nearly identical over 
the first 12 Gyr and only differ markedly in the most re- 
cent 1 Gyr. However deeper photometric views of lower 
luminosity Local Group dwarf galaxies find subtle differ- 
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ences between the morphological classes (Hidalgo 2011). 
The contrasting radial star formation histories between 
dlrrs and dSphs in that study may be due to differences 
in internal processes or environment dependent feedback, 
however disentangling the two and understanding why 
dlrrs retain more gas to power current SF is difficult. 

The presence of radial abundance gradients and their 
ubiquity within different classes of dwarf galaxies is one 
of the primary testing points to infer whether internal or 
external processes are responsible for the current abun- 
dance properties in dwarf galaxies. Depending on the 
mechanism for creating, sustaining, and erasing such 
chemical gradients, their presence and strength may cor- 
relate with physical properties o r environment of the host 
system (cf., Koleva et al.||20lT ). At present it is unclear 
if radial metallicity gradients are ubiquitous in low mass 
dwarf galaxies of all types, as seen in some dSphs flSa- | 
viane et al.|200lj|Harbeck et al.|2001||Tolstoy et al.|2004[ 



Battaglia et al.||2006[ |Kirby et al.||26l0| |Battaglia et al.| 
201 ip . Therefore it is also unclear if angular momentum 
and /or radial migration mechanisms mediate these gra- 
dients with different efficiencies in galaxi es with different 



dynamical histories, suc h as the dlrrs dSchroyen et al. 
2011] |Roskar et al.|[2008] ) 



There have been several large scale spectroscopic sur- 
veys of the RGB stars in the nearby [< 250 kpc) dSphs of 



the Local Group (e.g.jBa ttaglia et al 2006 Walker et al 



2009b 



Kirby et al. 2010). The gas rich dlrrs, lie at dis 



tances 500 — 1100 kpc away from the MW, which renders 
analogous surveys observationally expensive. Studying 
these isolated dlrrs is informative as they 1) provide im- 
portant initial structural and kinematic conditions for 
tidal trans formation scenarios which evolve dlrrs into 
dSphs (cf., |Kazantzidis et ah 2011a), and 2) being iso- 



lated offer a unique oppourtumty to study internal secu 
lar evolution in low mass halos that haven't been strongly 
perturbed by the tidal forces of the Milky Way. 
In the earli er papers in this seri es ( Leaman et al.|20"09 



Paper I) and ( Leaman et al.|2012 Paper 11) we presentee 



some of the first kinematic analyses of the stellar popu 
lations in an extremely isolated dlrr, WLM. This dwarf 
galaxy sits 1 Mpc from both the MW and M31, and 250 
kpc from the nearest neighbour, the low mass dSph, Ce- 
tus. From its Local Group position and velocity it is 
inferred that WLM's last pericentre passage was 11-17 
Gyr ago, which means it has h ad at most one close inter- 
action with a massive galaxy ( Leaman et al.|[2~012j ). The 
low mass and metallicity of WLM also make it an excel- 
lent laboratory for probing chemical evolution in regimes 



where the SF efficiency is expected to be low ( Krumholz 



& Dekel 2011 Kuhlen et al. 2011). WLM clearly has 



spent the majority of its lifetime m isolation, and there- 
fore is ideal for disentangling environmental and internal 
evolutionary processes. 

In this work we present the spectroscopic analysis of 
the Calcium II triplet (CaT) based [Fe/H] measurements 
in a sample of 126 red giant branch (RGB) stars in the 
dlrr galaxy WLM. As this represents the first isolated 
Local Group dlrr with a sizable spectroscopic survey of 
its evolved stellar populations, we proceed with compar- 
ing the dynamical, chemical, and structural properties 
to the less isolated dSphs and Magellanic Clouds. These 
comparisons are used to examine the interplay between 
environmental and internal feedback in an observational 



sense, and by selecting datasets that have substantially 
larger spatial extents and unbiased metallicity estimators 
than past studies, offer a significant improvement in the 
ability to accurately differentiate the chemical properties 
of Local Group dwarf galaxies. 

2. OBSERVATIONS AND DATA SOURCES 

Resolved stellar spectroscopic data for this paper are 
from the following sources: spectra of 180 RGB stars 
were observed with Fo cal Reducer and low dis persion 
Spectrograph (FORS2, Appenzeller et al. 1998[ ) at the 
VLT (Paper I) and t he Deep Imaging Multi-Object Spec- 
trograph (DEIMOS, iFaber et al.|2"003| at Keck II (Paper 
II) for WLM. Metallicities from FORS2 data were pub- 
lished in paper I, and velocities from both samples were 
published in paper II. In this paper we determine metal- 
licities from DEIMOS data, as well as re-calibrate the 
FORS2 EW obtained in paper I. The dSph RGB spectro- 
scopic data come from the following sources: the Sculp- 
tor, Fornax, Sextans, and Carina data are taken from the 
original surveys of 



( [20061 [20TT| >; [Koch 




pies have since been updated by 
and Radial velocity Team (DART). The updates include: 
observations of additional stars with the same instru- 
ment setup and reduction as described in the original 
paper J4 applying the new CaT- [Fe/H] calibration from 
Starkenburg et al. (20101 to all stars, visual inspection 
of low metallicity candidates, and quality control cuts 
(S/N> 10, Svhei < 5 km s _1 ). Throughout this paper 
spectra from these four galaxies will be referred to as 
the "DART sample" , meaning these updated catalogues 
which are based on the original papers. The full sample 
of eight dSphs from [Kirby et al. | ( |2010[ ) (hereafter K10) 
were also considered, however as the goal was to char- 
acterize the chemical properties of the complete spatial 
extent in these dwarf galaxies, only Leo I and Leo II are 
used from that work, as they both show spatial coverage 
out to at least three core radii (and past the tidal radii). 

The WLM and dSph data are supplemented with ob- 
servations of RGB field stars in the Magellanic Clouds . 
For the LMC we draw f rom the work of [ Cole et al.| ( |2005[ ) ; 
Pompeia et al.| J2008(); iCarrera et al.| ([2lio5rj|anT The 
studies of Carrera et al. ( 2008a )~ |Parisi et ai.| ( |2010[ ) for 
the SMC. The data taken from literature represent spec- 
troscopic surveys of at least 100 stars per galaxy where 
[Fe/H] estimates are available for each star, and the sam- 
ple extends to at least 75% of the tidal radius, r t , of the 
galaxy. In addition WLM and most of the DGs compiled 
in this work are sampled well along both the minor and 
major axes, providing a high degree of angular unifor- 
mity. The data sources and references for the nine Local 
Group dwarf galaxies used in this paper are summarized 
in Table 1. We refer the reader to the individual stud- 
ies for more detailed information on the reductions, but 
throughout this paper differences in methodology will be 
discussed when relevant to comparisons we draw. 

2.1. Auxiliary Measurements 



1 The exception is Carina, whose stars from [Koch et ah[(|2006| 
were taken from the ESO archive and reprocessed using the IJArtl' 
pipeline procedures. 
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Additional parameters aside from the published, di- 
rectly measured velocities and metallicities are required 
for our analysis. These are core radii of the dwarf galaxy 
(r c ), elliptical radii and ages of the stars. With these it 
is possible to explore the spatial and temporal variations 
that may illuminate differences or connections between 
the evolutionary history of dlrrs and dSphs. For the 
elliptical radii measurements, a galaxy centre, mean el- 
lipticity (e) and position angle (PA) are required. In the 



case of WLM the field center from Mateo ( 1998 ) , and 
the ellipticity, position angle, and r c from Paper 11 were 
adopted. For Leo I and Leo II dSphs we adopt the fi eld 
centre from K10 and positio n angle, r c and e from |Ir-| 



win & Hatzidimitriou ( 1995 ) .Values for the field centre 
for the DART dSphs were taken from the original papers 
listed in Table 1 and references therein, and core radii, 
e, and PA from 



centres, e and 



Irwin fc Hatzidimitriou] ( |1995 ). T he field 
J A lor the SM C were taken fr om Harris | 



fc Zaritsky l ( |2006[ ), and from |van der Marel| ( |2001[ ) for 
the LMC. The r c for the LMC and SMC were computed 



mputed 

as a weighted mean of the values in the recent study 
of Belcheva et al. (20111. The geometrical radii for all 
samples have been derived, and in each case these have 
been referred to the global ellipticities and position an- 
gles listed above. 

3. WLM SPECTRAL ANALYSIS 

For the spectral analysis of WLM we have used the pre- 
viously observed FORS2 spectra as well as the newly ac- 
quired DEIMOS data. In merging the two datasets pos- 
sible systematic offsets due to instrumental signatures, 
equivalent width (EW) measuring techniques and metal- 
licity calibration must be explored. The following section 
outlines the influence of the metallicity calibrations on 
the full sample, with special emphasis on joining the old 
and new datasets in a way that is consistent within the 
errors. This allows for a homogeneous joint sample to be 
created which has global metallicity properties that are 
insensitive to instrumental and calibration details. 

As in Paper I, [Fe/H] values were derived based on the 
empir ically calibrated Calcium II tripl e t (CaT) method 
'e.g. 



Armandroff & Da Costa 



19971. Originally calibrated with old galactic globular 



1991 Rutledge et al. 



clusters, work has also bee n done extending the CaT 



calibration to younger ages (Cole et al. 20041, which is 
relevant here given that the R UB popu lations in WLM 
are expected to span 10 Gyrs (| Dolphin]|2000 [). However 
due to the fact that the empirical calibrators (Galactic 
globular clusters) are only found as metal poor as [Fe/H] 
~ —2.3, the method has intrinsic limitations. In recent 



years work by Kirby et al. (2008) has shown the advan- 
tage of using synthetic spectral techniques that can de- 
rive [Fe/H] in a way that is not limited by the properties 
of the calibrators - thus pushing down to lower metal- 
licities. Additionally, the empirical li near calibration of 
the Ca T method has been revised by |Starkcnbur g~et al.| 
(2010) in order to address the CaT calibration limitation. 
They found that at low equivalent widths and magni- 
tudes, the CaT behaviour becomes non-linear, resulting 
in an overestimate of the [Fe/H] for low CaT equivalent 
widths. 

3.1. Equivalent Width Measurements of DEIMOS 
Spectra 



Equivalent width measurements of the Calcium II 
triplet lines in the new DEIMOS spectra were done with 
pixel-pixel integration methods, as in Paper I. For com- 
parison, equivalent width measurements were also pro- 
duced using integration over the line and continuum 
bandpasses of |Cenarro et al. (2001 1, and fits to the lines 
using Gaussian, Lorenztian, Moffat and Voigt functions. 
A comparison of the pixel integration versus several fitted 
equivalent width estimates for stars in the DEIMOS cali- 
bration clusters shows exc ellent consistency, ho wever the 
integration bandpasses of Cenarro et al. (2001) tends to 
produce larger values for the KW relative to other meth- 
ods. For comparison, the difference in EWs measured 
by pixel integration and Gaussian plus Lorentzian fits 
in the lower resolution spectra from Paper I, translated 
into a 0.17 dex difference in [Fe/H] - which as we will 
demonstrate is much less than our random uncertainties. 
For these data which have relatively low spectral resolu- 
tions the direct pixel integration shows the best agree- 
ment with Gaussian fits for the calibration stars over 
all metallicities. However as in Paper I, the resolution 
and signal-to-noise of the WLM stars necessitate using 
the integrated equivalent widths rather than functional 
fits, as the line FWHMs are on the order of the spectral 
resolution and contaminating noise features nullify any 
difference between the functional fits. 



3.2. Placement of the Joint Sample onto the Metallicity 

Scale 

We derive metallicities for WLM using five different 
CaT- [Fe/H] calibrations to explore th e variation in globa l 
metallicity prope r ties, t hese include: Cole et aL] (|2004[) 
i Battaglia et al. (2008), two from Starken burg et ah] 
( |2010| , and one based on our DEIMOS calibrating clus- 
ters. Any of these CaT- [Fe/H] calibrations requires a 
summed equivalent width determination for each star, 
and are most precise when the EW measurements are 
done using the same line measur i ng tec hnique as the orig- 



inal calibrations. In Cole et al. (2004), each of the three 



calcium triplet line measurements were combined in an 
unweighted fashion to yield a summed equivalent width 
per star YW = W si9& + W S5 4 2 + W 8e62 - With this re- 
lation the calcium index W = YW + ft(V — Vhb) is 
formed. The term in the parentheses provides a cor- 
rection for the changes in T e ff and log(g) for stars in 
different phases on the red giant branch. Our V mag- 
nitudes for WLM are taken from the INT WFC cata- 



logue presented in |McConnachie et al. ( 2005 ) and the 
calibrations in that paper. We adopte d the horizon- 
tal branch at Vhb = 25.71 ± 0.09 mag (Rej kuba et al. 
0.73±0.04A mag 



2000), and take 



from 



Cole et al 



fl2004[ ). T he photometric studie s of |McConnachie et al 
2005p and Rejkuba et al. (2000) find distance moduli ot 
m - M) = 24.85 ± 0.08 an7T(m - M) = 24.95 ± 0.13 
respectively, indicating that there are no significant zero- 
point offsets that could influence our CaT- [Fe /H] calcu- 
lations. Using the Carretta & Gratton (1997) scale, the 
calcium index is converted to a metallicity (]Fe/H]cc797) 
using Equation 3 of Paper I. 

The summed equivalent widths were al so converted to 
an [Fe /H] scale using the calibration from Battaglia et al. 
([20 08]), as well as two non-li near calibrations presented 
in ( jStarkenburg et al. 20101. In these three cases the 
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summed equivalent width was computed using the two 
longer wavelength calcium triplet lines. We also created 
a linear calibration based on our stellar cluster (GC) cal- 
ibrators (NGC 6791, Pal 14, NGC 7078) to check that 
there were no strong dependencies on instrument resolu- 
tion. 

Figure [T] shows the summed equivalent width as a func- 
tion of V magnitude relative to the Horizonta l Branch 
for the calibration based on Equation 5 from Starken- 
burg et al.| Jg Olo]). The summed EWs of the DEIMOS 



stars track the Starkenburg et al. (2010) relations closely, 
perhaps suggesting that at least for moderate resolution 
spectra, pixel int egration of the CaT lines produces EWs 
which fall on the Starkenburg et al. ( 2010 1 relations. The 
non-linearity in the low EW and taint end of the param- 
eter space is evident in this diagram. While the calibrat- 
ing Globular Cluster stars shown in that figure extend 
to low magnitudes below the horizontal branch, we note 
that most of our WLM member stars are in the region 
of the parameter space, where large differences between 
the linear and non-linear calibrations are not expected. 

While there is a more recent [Fe/H] sca le based on glob- 
ular clusters from Carretta et al. ( 2009[), it is an average 
of fou r past scales (including the (Jarretta & Gratton 
(1997) scale). Unfortunately only one of our calibrat 



( 2009| sam ple (NGC 7078) for which the |Starkenburg 



et al 



ing clu sters is directly measured in the JCarretta et al 



the F< 



= p|2010[ ) calibration shows excellent agreement 
' ORS2 calibr ating clusters alon e we find that 1 



From 
the dif- 



ference using the Cole et al. (|2004) calibration between 



the Carretta & Gratton 1 TW7) ilmcflCarretta et al. (12009 1 
scales for JNGC 104 and NUU 7078 is only 0.06 aAcTO^D 
dex respectively. Given the relative size of the random 
uncertainties, and other factors discussed in joining these 
merged datasets this choice of absolute [Fe/H] scales does 
not introduce a change in our analysis. 

This is apparent from examining the metallicity dis- 
tribution functions (MDFs) for 126 stars of the 180 
member stars which had sufficiently high signal-to-noise 
(> 10A _1 ) - shown for each of the five calibrations in 
Figure [2] While the calibrations are all in good agree- 
ment with one another for the relatively bright and metal 
rich stars of WLM, we wish to compare WLM to low 
metallicity dSph sy stems and therefore adopt the non- 
linear calibration of |Starkenburg et aT| ( |2010[ ) based on 
the horizontal branch magnitude tor the analysis in the 
rest of this paper. This allows us to compare WLM to 
the dSphs on a consistent scale that does not suffer from 
saturation biases which would be prevalent in the faint, 
low metallicity dSph stars. 

Metallicity calibration biases are especially important 
to consider when studying the spatial distribution of 
metallicities within a galaxy. In Figure [3] we show the 
effect on spatial variations in [Fe/H] for two different 
CaT based [Fe/H] calibrations for WLM. There are vis- 
ible differences even for the relatively bright and metal 
rich stars in our sample, which would be even more se- 
vere in low metallicity systems. As gradients have such 
large implications for dwarf galaxy formation and evolu- 
tion, care must be taken during interpretation, and this 
starts with an accurate understanding of the calibration 
biases. 

Uncertainties on the metallicities were propagated 



from the initial line width measurements to the cali- 
brated [Fe/H] values, as in Paper I. The mean uncer- 
tainty of the FORS2 stars is 0.25 dex, and 0.26 dex for 
the DEIMOS stars. While higher resolution, some of the 
DEIMOS data are much lower signal to noise, which al- 
lowed us to derive reliable metallicity estimates for only 
50% of those spectra. 

3.2.1. Consistency Checks 

The MDFs presented in Figure[2]show good qualitative 
agreement and have similar dispersions and mean values 
within the uncertainties. While ideally a homogeneous 
instrument setup for the full sample of stars is preferred, 
we find consistent metallicity signatures between Paper 
I and this joint sample within the large uncertainties. 
Evidence for this comes from repeat measurements 
of the two stars in common between the DEIMOS 
and FORS2 observations (which show a difference of 
0.13 < A[Fe/H] < 0.31), and a similar mean metallicity 
(-1.27 ±0.04; -1.28 ±0.02) for all of the stars in Paper 
I and this full sample. Additionally the radial gradient 
computed in Paper I (d[Fe/H]/dr r ~ -0.076 ± 0.03 



dex r 1 ) is unchanged when we apply the Cole et al. 
(2004) calibration to the joint sample of FORS2 and 
DEIMOS spectra (d[Fe/H]/dr c ~ -0.08 ± 0.03 dex 
r^ 1 ). Therefore, while there are subjective choices 
on the metallicity scale, EW measurements, and CaT 
calibration, these changes are all within the errors. This 
suggests that given the large uncertainties inherent in 
spectroscopy of stars in distant galaxies like WLM, our 
joint MDF is sufficiently robust to analyze the global 
metallicity properties 



3.3. Age Derivations 

Age derivations were discussed in Paper I and Paper 
II for the WLM sample. Ages were derived using the 



published photometry and the Demarque et al. (20041 



stellar evolution models. The older library is chosen due 
to the metallicities in the larger sample and new cali- 
bration being outside the range of the Victoria-Regina 
models used in Paper I. In addition a greater flexibility 
in a— abundances is possible with the Demarque et al. 
(2004) models. The V and I photometry, reddening, 
and distance moduli were taken from Papers I and II on 
WLM, and discussed therein. The ages were interpolated 
using the grid of isochrones, the dereddened photometry, 



and spectroscopic 



Paper II the 
tion of [Fe/H 
eTaD pOlOl; 



a/Fe 



Fe/H] and [a/Fe] estimates. As in 
values were interpolated as a func- 
using the l i teratu r e values from Colucci 
Venn et al.| ( [2003] ); |Bresolin et al | ( j2006| ) 



to describe the mean trend of [a/Fe J versus [Fe/HJ in 
WLM. Errors were assigned by propagating the photo- 
metric, reddening and distance modulus uncertainties, as 
well as the spectroscopic abundance uncertainties into 
the position of the star on the colour magnitude dia- 
gram. Ages derived using this method will be valid in 
a differential sense within a sample, as there are strong 
systematic uncertainties between the stellar evolution li- 
braries used in various studies. However, the metallicity 
uncertainties dominate over the choice of evolution li- 
brary for such distant systems as WLM, therefore the 
general age-metallicity properties of WLM may be ex- 
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tracted. The relative random uncertainty on age for an 
individual WLM star is ~ 50%. 

3.3.1. Quantifying Systematic Age Errors 

The random error captures the uncertainty in derived 
age due to errors in colour, magnitude and [Fe/H], how- 
ever there are three additional systematic errors not in- 
cluded in the previous section that must be quantified 
- AGB contamination, differential (internal) reddening, 
and v ariations in [q /Fe], WLM exhibits an extended 
SFH ( |Dolphin|[2000l ), therefore it is highly probable to 
sample stars on the giant branch with ages 1.6 — 12 Gyr. 
In addition the distance of WLM makes it difficult to 
accurately differentiate second ascent giant branch stars 
with much confidence from photometry. This means that 
within the sample there may be AGB stars; these do not 
affect the derived [Fe/H] or velocities but can produce 
a bias in the inferred age. Using the SFH of [Dolphin 
( |2000[ ), it is possible to roughly estimate the AGB con- 
tamination rate within the colour and magnitude range 
of the WLM spectroscopic targets. A conservative upper 
limit on the contamination fraction is ~ 40%, with about 
1/3 of those AGB stars being younger than 2.5 Gyr, and 
a third older than 9 Gyr. Using a grid of isochrones, it is 
possible to work out for a given colour and magnitude the 
difference in age between an RGB and AGB star. The 
systematic age error due to AGB contamination is found 
to be strongest at young ages. Specifically, an AGB star 
of 1.6 Gyr would have its age underestimated by 20% 
if it were considered an RGB star in the sample. This 
percentage drops to 10% for a 2 Gyr star, and 5% for a 
10 Gyr star. 

The unknown nature of differential internal reddening, 
and star-to-star variations in a— element abundances in 
WLM stars will contribute additional systematic errors. 
To numerically estimate the combined systematic uncer- 
tainty due to these two factors and the above mentioned 
AGB contamination, we proceeded as follows. For a test 
star of a given true age and [Fe/H], and fixed evolution- 
ary position ^0.5 magnitude below the tip of the RGB, 
we drew randomly a possible variation in (V-I) (due to 
reddening), and [a/Fe], as well as gave it a 50% chance 
at being an AGB star. The distribution of internal red- 
denings was drawn from a Gaussian of ay-i = 0.03 mag 
(for comparison, the line of si ght reddening in the direc- 
tion of WLM is E(V-I)= 0.037jSchlegel^t^ni998;). The 
distribution of differential [a/Fe] was taken from a Gaus- 
sian of a a — 0.05 dex - which was chosen primarily to 
keep [a/Fe] within the range of the isochrone grid. In a 
given iteration, the test star had its age rederived using 
the new colour and magnitude on a grid of isochrones re- 
flecting its new [a/Fe]. If the star was also drawn to be 
an AGB star, the new age was modified by the systematic 
age offsets discussed in the above paragraph. 

In Figure H] we show the systematic errors due to these 
combined effects for 10000 iterations on each input star 
of a given canonical age and [Fe/H]. While AGB con- 
tamination most strongly impacts the age systematics 
at young ages, the effect of differential reddening and 
a— element variation dominates for low metallicity stars. 
The standard deviation of Monte Carlo trials are indi- 
cated as ellipses, with the mean movement indicated by 
the black arrows. Where appropriate in this work and in 
Paper II, we adopt the semi-major axis of the ellipses as 



an estimate of the total systematic uncertainty for stars 
of various ages and metallicities, and incorporate that 
along with the individual random uncertainty on age for 
a star. The distance of WLM makes deriving ages dif- 
ficult, however we are aided by the relatively metal rich 
and young populations of stars in this sample. 

4. RESULTS 

We now compare the distribution of metallicities, both 
spatially and with age (where possible), using the spa- 
tially extended sample of metallicities for hundreds of 
stars in each of WLM, LMC, SMC, Fornax, Leo I, Sculp- 
tor, Leo II, Sextans, and Carina. Through this analysis 
we will test whether an isolated galaxy such as WLM, 
with its relatively quiescent tidal evolution and gas con- 
tent history, shows differences in chemical evolution com- 
pared to the gas poor-tidally disturbed dSphs or gas rich 
tidally disturbed Magellanic Clouds^ 

4.1. Metallicity Distribution Functions 

In comparing the metallicity distribution functions 
(MDFs) between the dwarf galaxies, it is important to 
be aware of the sample sizes and spatial extents of the 
datasets for comparison galaxies, and where [Fe/H] has 
been calculated using differing methods. Both the K10 
and the DART sample have [Fe/H] estimates that should 
be free of calibration biases (see §3). There remains sys- 
tematic differences due to the varying spatial coverage 
and number of stars in each study however. The global 
metallicity distribution may be biased in studies where 
only central regions of a galaxy are sampled, especially 
in the case where there is a radial metallicity gradient. 
Similarly when the sample size is small the metallicity 
may be underestimated, as a minority metal poor pop- 
ulation is more difficult to sample as efficiently due to 
a population bias resulting from mass dependent stel- 
lar evolutionary timescale^j This effect will be enhanced 
in the presence of any age-mctallicity relationship and 
is still present for populations with extended SFHs like 
WLM. 

To study the impact of these biases, MDFs were com- 
puted for the galaxies using an equal spatial range, and 
with equal numbers of stars. The samples were restricted 
in these test cases to stars within 1.5r c of the dwarf 
galaxy center (roughly the smallest radial extent of dSph 
samples in the K10 catalogue). Any one galaxy that has 
a strong metallicity gradient will be biased in this case, 
however the MDFs of the inner regions will be much more 
appropriate to compare to each other in a differential 
sense. To study the bias in sample size, the datasets 
were also resampled from the inner 1.5r c so that they 
have equal numbers of stars. In this case the sample 
sizes have been constructed to have 31 stars (equal to the 
smallest sample of stars within 1.5r c in any of the galax- 
ies) . Each of the galaxies were resampled from the larger 

2 It should be kept in mind that tidal influences on the Magel- 
lanic Clouds may have on ly begun recently (< 2 Gyr) if they are 



on fir st infall to the MW dBesla et al. 12007 Boylan-Kolchin et all 
[20111 1 ' iLJL ' 

3 For a given time interval, more young, high mass, high metal- 
licity stars will evolve to the TRGB than older, low mass, metal 
poor stars in the same time interval. The dominant fraction of 
young stars at the TRGB is large enough to outnumber the metal 
poor old star s, despite the con tradictory initial relative numbers of 
the IMF. cf., |Cole et al.|2008| 
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population 10000 times, computing a cumulative metal- 
licity distribution function (CDF) each time. An average 
CDF was constructed based on 31 stars in each galaxy. 
Figure |5] shows an example of the mean (solid blue) and 
dispersion (dotted blue line) of rcsampled CDFs for sev- 
eral of the galaxies in this study, as well as the impact 
of the radial cutoff. This shows that the primary bias 
affecting comparison of different studies is their spatial 
extent, rather than any sampling bias from differing num- 
bers of stars. We have therefore selected surveys which 
have resolved spectroscopic data out to at least 75% of 
the galaxy's tidal radius in hopes of minimizing any spa- 
tial bias when comparing the dwarf galaxies, as shown in 
Table 1. 

Figure [6| shows the MDF of WLM, the Magellanic 
Clouds and the six dSphs, ordered by luminosity. The 
full sample and central regions are shown for each galaxy. 
In each panel the 10th, 50th, and 90th percentile metal- 
licities of the MDFs are also shown (as dotted lines). 
The mean (median) metallicity of WLM is [Fe/H] = 
— 1.28(— 1.24) ± 0.02 dex, in good agreement with the 
trend shown by the Local Gro up luminosity-metallicity 
relation (cf., Woo et al. 2008). The uncertainty on the 
value for WLM represents the error on the mean, however 
typical systematic uncertainties for both synthesis and 
equivalent width based measurements from the literature 
are ~ 0.15 dex. The median value of WLM lies within 
the range of [Fe/H] distributions shown by the dSphs, 
and is in close agreement to the more luminous members, 
Fornax and Leo I as well as the SMC. The metal poor 
population below [Fe/H] = —1.74 forms a ~ 10% minority 
population in WLM, and stars are found as metal poor 
as [Fe/H]= —2.85 dex. RGB stars with metallicities as 
enriched as [Fe/Hl= —0.35 dex are found, consistent with 



the results from Venn et al. (12003) for two supergiants 
in WLM. Table 2 tabulates the percentiles of the MDF 
for each galaxy. The mean and extrema metallicities of 
the isolated, gas rich WLM dlrr show little differences 
from the luminous dSphs, nor the gas rich but (recently) 
tidally perturbed SMC. 

WLM shows a similar median metallicity to some clas- 
sical dSphs and the SMC, however there appears to be a 
slight preference for the lower mass dSphs to show more 
asymmetric, extended metal poor tails relative to their 
mean metallicity. Fornax, Leo II, Sextans, and Carina all 
show a significant asymmetry to their MDF compared to 
the dlrrs - however it may be that the metal poor stars 
in the gas rich dwarfs are simply more difficult to sample 
efficiently in such cases with ongoing star formation (see 
§4.1). The difference in metal poor tails can be seen qual- 
itatively in the right panels of Figure [6j where leaky box 
and pre-enriched chemical evolution model s have been 
overlaid. The formalism of |Prantzos| ( |2008[ ) for the two 
models has been used, and the best fitting models deter- 
mined through a maximum likelihood approach assuming 
Poisson errors on the distributions. The estimated yields 
and initial metallicity for both chemical evolution mod- 
els are listed for each galaxy in Table 2. Qualitatively 
the simple leaky box models provide reasonable fits to 
many of the dSphs, but the dlrrs and Leo I appear more 



the DART sample of Fornax (which is nearly the same 
luminosity as Leo I) samples out twice as far in radius 
and s hows a significant metal poor tai l. Alternatively for 
Leo I, Lanfranchi & Matteucci (2010) showed that mod- 
els with intaliing pristine gas could also reproduce the 
MDF. 

The spread in [Fe/H] of a dwarf galaxy has been found 
to anticorrelate with mean metallicity in the recent study 
of Kirby et al. (2011 ). We revisit this relation here, as our 
literature sample explores both higher luminosity galax- 
ies (some of which are gas rich, and one of which is iso- 
lated), and typically have larger spatial coverage. The 
intrinsic dispersions have been calculated by subtracting 
in quadrature the mean error in metallicity for a sample 
from the total measured dispersion. Figure [7] shows the 
intrinsic spread in metallicity, (J[ Fe/H]_i f° r the aata in this 



paper as well as those derived in Kirby et ahl ((201 1|) and 

TT 



the linear correlation those authors found. The uncer- 
tainty on th e dispersions are calculat ed using relations 
from §3.1 of Hargreaves et al. (1994) and the intrinsic 
and raw dispersions from this work and K10. Figure [7] 
suggests that the K10 relation still provides a good de- 
scription of the large spread seen in [Fc/H] for the Ultra 
Faint Dwarfs (altho ugh Bootes I wa s recently found to 
lie off this relation; Lai et al. 2011). However at high 
luminosities (> 1O 5 L ) t he dispersion in me tallicity may 
saturate, as hinted at in 
little change even up to 



Norris et al. (2010), as it shows 
he brighter dlrrs such as WLM 
and the Magellanic Clouds. This may be due to the high 
luminosity systems having an enrichment timescale that 
is similar, or above some thres hold duration t o produce 
a well mixed ISM (see §5.1; also |Leaman|20 1 2 ) . Alterna- 
tively it could suggest that the lowest mass systems are 
more impacted by gas and metal outflows, driving them 
to larger metallicity dispersions, compared to the higher 
mass galaxies which retain a larger fraction of their gas 
and metals. 



4.2. Spatial Variations in Chemistry 

Paper I showed that the metal rich and metal poor 
stars in WLM are similar in their spatial distributions 
with slightly more metal rich stars in the inner regions. 
With only 78 stars in the original WLM sample how- 
ever, it was difficult to interpret the change as due to a 
superposition of two populations, or a smooth gradient. 
Here we analyze the spatial distribution of the WLM stel- 
lar populations using the new DEIMOS data and CaT- 
[Fe/H] calibration, and examine the metallicity gradient 
in comparison to tho se of dSphs (i.e., Tolstoy et al. 12004"* 



|Ba ttaglia et al.|2006 ) . With the new stars and calibration 



closely fit by the p re-enriched solutions (see also |Gul- 
lieuszik et al. 2009). In the case of Leo I this could be 
due to the low metallicity stars lying at large radii (see 
§4.2) outside the spatial coverage of the K10 survey - as 



by |Starkenburg et al. (2010), the metallicity distribution 
has changed and the spatial signatures are altered from 
Paper I. Figure [3] shows that the new calibration clearly 
favours a milder gradient. 

Metallicity gradients in d Sphs are sti l l com mon to 
varying degrees as shown in Kirby et al. (2011). How- 
ever the spatial coverage in that work is relatively lim - 
ited compared to the samples from Walker et al. (2009a), 
or the DART survey. Here we can explore compar- 
isons of [Fe/H] gradients in a rigorous manner for the 
Local Group galaxies due to the large spatial coverage 
{r max > 0.75r t ) in the literature sample. Figure [8] shows 
[Fe/H] as a function of elliptical radius (in units of core 
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radii), with the dashed lines showing linear least squares 
weighted fits to the data. 

WLM shows a mild gradient with d[Fe/H]/dr c = 
—0.04 ± 0.04 dex r" 1 , similar to the Magellanic Clouds. 
It should be noted that the LMC RGB metallicities of 



Carrera et al. (2008b I show a correlation with magni- 
tude, which manifests itself as a systematic steepening 
of the gradient. Correction of this reduces the measured 
gradient by half, consistent with an extended sample of 
RGB stars from a forthcoming study (d[Fe/H]/dr c = 
-0.008 ± 0.004 dex r" 1 ; Cole et al., in prep.) - there- 
fore the LMC gradient shown here could be taken as a 
lower limit. In contrast the dSphs show radial metallicity 
profiles which are much steeper. 

Many of the individual galaxies in Figure [8] show 
complex trends of [Fe /H] with radius, and as noted in 
Battaglia et~aL"1 ( |2011| ) some galaxies show steep profiles 
in the inner few core radii, which then flatten to a low 
metallicity plateau. In such cases the gradient may be a 
superposition of two populations with differing concen- 
trations - a possibility that is difficult to rule out quanti- 
tatively. We have computed running boxcar averages of 
[Fe/H] versus radius in Figure[9]to track the mean metal- 
licity more precisely. The dlrrs still show statistically 
flatter radial [Fe/H] gradients relative to the dSphs (es- 
pecially considering the uncorrected systematic with the 
LMC profile). The low metallicity plateau in Sculptor 
is quite evident at large radii, and while there is scatter 
between any given dSph, overall they show much steeper 
metallicity drop offs in the inner three core radii. Fig- 
ure [9] clearly shows the necessity for very spatially ex- 
tended spectroscopic surveys in order calculate the com- 
plete abundance gradient and global metallicity distri- 
bution, as within the inner 1.5r c of the dSphs there is 
a stochastic behaviour before the abundances coherently 
fall at larger radii. 

4.3. Age Metallicity Relations 

Photometric SF histories already suggest that WLM 
experienced ex tended star formation over its lifetime 
(Dolphin 2000), and with spectroscopic [Fe/H] and age 
estimates in WLM, we can compare the age metallicity 
relation to several othe r Local Group dwarfs in a differ- 
ential sense. Figure 10 plots the age metallicity relations 
(AMRs) for the three dwarf irregulars as well as For- 
nax. For the LMC and SM C we have used the publ ished 
age-metallicity rela tions in Cole et al. (2004) and Car- 
rera et al. (2008a I respectively. Similarly with Fornax 
we have ad opted the published m ean age-metallicity re- 
lation fromlBattagliaeFal. (2006). In those three studies, 
stars with ages older than 10 Gyrs have been scaled as 
tnew = 10 + 0.41(f — 10), in order to homogenize the 
oldest stars to a consistent maximum age, as the studies 
differ in the uppermost age the stellar evolution libraries 
consider. For WLM five bins have been computed for 
the AMR, and the error on the mean for the metallicity, 
and combined random and systematic age errors com- 
puted for each bin (§3.3.1). Where applicable the field 
RGB star data are supplemente d with the [Fe/H] val - 



ues from supergia nt s tudies by |Bresolin et al.| ([20061) 

imilarly 



Venn et al 
the oldest 



(2003) and Levesque et al. 1(2006 



WLM flColucci et al.(2010| ) were added and metallicities 
taken from those spectroscopic studies. There appears 
to be good agreement between these anchor points and 
the youngest and oldest field RGB star age-metallicity 
data. The simple leaky-box and pre-enriched chemical 
evolution models which were fit to the MDFs in Figure 
[6] have been overlaid as well. 

The AMRs of SMC, Fornax, and WLM in Figure M 
are all qualitatively similar in their shape, with relatively 
shallow metallicity enrichment over intermediate ages, 
and only an offset in metallicity at each age. WLM shows 
slightly more rapid metal enrichment in the last 3 Gyr 
compared to any of the other galaxies, consis tent with the 
burst of central SF found by Dolphin ( 2000[ ). The nearly 
flat [Fe/H] values from 3 — 9 Gyrs in WLM are quite 
similar to the SMC for that time frame, but Fornax shows 
a stronger increase in metallicity over the same timescale. 
This is likely due to the strong burst of SF estimated to 
have occurred bet ween 3-8 Gyr ago in Fornax based on 
its CMP analysis ([Stetson et alLl998 ' 



Tolstoy et"aLl|20011 |Coleman & de 



Saviane et aT[2 000 
ong||2008t Ide Boer 



|et al.||2012p . For the first two Gyr of it's lifetime, WL 
does not appear to have undergone as rapid enrichment 
as the LMC. The mean metallicity of WLM changes by 
< 0.5 dex in this period, similar to the SMC and Fornax, 
while the LMC changes its mean metallicity by nearly 1 
dex. 

The AMR of the LMC and SMC has be en carefully 
examined by |Pagel &: Tautvaisiene ( 1998 1 , where the 
primary difference in their AMR shapes was attributed 
to the relative strength of the SF burst (s) in the LMC. 
WLM's SFH has be en interp r eted fr om HST colour mag- 
nitude analysis by Dolphin (20001, and the computed 
AMR from that work is also shown overlaid on the WLM 
data in Figure [To} There is remarkable agreement be- 
tween the shape of the AMR derived from our spectro- 
scopic s ample, and that inferred from the photometric 
study of iDolphinl (|2000b. The similarities in the AMR 



of WLM to that or the dSph Fornax, illustrate that an 
isolated dlrr like WLM (which has experienced much 
less tidal perturbations over its lifetime, and exhibits dif- 
ferent kinematics) has had chemical enrichment proceed 
similarly to both gas rich but tidally perturbed (SMC), 
and gas poor, low angular momentum (Fornax) systems 
of similar luminosity. 

The dispersion in metallicity for stars of differing ages 
can provide insight into such processes as radial migra- 
tion. For the youngest stars in WLM (< 2 Gyr) the 
intrinsic dispersion in [Fe/H] is found to be a±—2 = 
0.15 ± 0.04 dex. This increases to cr 3 _ 9 = 0.32 ± 0.07 
dex for stars 3 — 9 Gyr, and <7io-i3 = 0.33 ± 0.09 dex for 
stars > 10 Gyr. Using the sample of Cole et al. (20051 



and identical age cuts for the LMC, the intrinsic disper 
sion in [Fe/H] is a^ 2 = 0.13 ± 0.01, <t 3 _ 9 = 0.25 ± 0.02, 
and crio-13 = 0.60 ± 0.10. For both the LMC and WLM 
there appears to be an increase in the metall icity spread 
for old er ages, however as pointed out by |Cole et aL 



obular clusters in the LMC (|Colucci et al 



2011| ), Fornax ( |Strader et aT||2003[ ) and the one GC in 



([20051 this can be due to systematic age effects blurring 
the age-metallicity relation. In the case of the LMC the 
rapidly rising [Fe/H] enrichment in the first few Gyrs also 
adds to the high intrinsic spread in [Fe/H] at large ages. 
However the AMR is relatively flat in early times for 
WLM, so the increased [Fe/H] dispersion at large ages 
could be linked to a radial mixing process. Both the 
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LMC and WLM also show an increase in velocity disper- 
sion with age, as well as increase in metallicity dispersion 
with age, however modelling of such joint chemodynamic 
signatures is out of the scope of this paper. 

5. DISCUSSION 

In the following section we asses the likelihood that 
secular or environmental factors could produce the par- 
ticular similarities and differences observed in the metal- 
licities, structure, and dynamics of the sample of galaxies. 

5.1. Global Metallicity Properties in the Sample 

Despite WLM's isolation, Figures [6] and [7] suggest that 
its bulk chemical properties are similar to the more lu- 
minous of the tidally disturbed dSphs and the Magel- 
lanic Clouds. Approximately the same mean metallicity 
is found in WLM, Fornax, Leo I and the SMC, while 
their tidal indice^j range from O = 0.2 to = 3.?^ 
( Karachentsev|[2005| ) . Stars are found as metal poor as 
[Fe/HJ~ —2.9 m WLM, the same lower limit as Fornax, 
suggesting that dSphs and dlrrs of the same luminosity 
have star s of low metallicity p resent - consistent with the 
survey of Weisz et al. (2011a). Together with the similar 
AMR shown lor these galaxies in Figure [lOj the com- 
parative metallicity properties suggest that the chemical 
enrichment proceeds largely independent of environment 
and is primarily dictated by the m ass of the galaxy 



agreement with the simulations of Sawala et al. (2011). 
This is also consistent with the recent observations of 



the isolated dwarf galaxy VV124 by jKirby et al.| ( |2012[ ), 
which was found to lie on the Local Croup MMK despite 
its extreme isolation. 

Similarities in the AMR and chemical enrichment 
timescales of the SMC, WLM and Fornax may also be 
indicative of a common initial mass - with subsequent 
differences in the present day masses plausibly attributed 
to tidal stripping by the M ilky Way. Using e stimates for 
the half light masses from Wolf et al. (20101 and Paper 
II, in conjunction with the tidal evolutionary tracks of 
Penarrubia et al. (2008), it can be shown that if the SMC 
were to have lost ~ 40% and ~ 98% of its total mass, it 
would have a similar half light mass as WLM and Fornax 
respectively. Similarly were WLM to undergo stripping 
of ~ 90% of its total mass, it would have a compara- 
ble present day mass to Fornax. It is certainly plausible 
then that the similar enrichment history of those three 
galaxies is linked to similar infall masses. Adding weight 
to this argument is the presence of GCs in both WLM 
and F ornax, as well as plan etary nebulae ( |Magrini et al. 



2005 Kniazev et al. 2007 [) , which may imply common 
masses and SFH (cf.' paviane et al.||2009| ). 

The measured global metallicity properties of WLM 
are in good agreement with the expectations for its lu- 
minosity, regardless of environment. WLM shows no dis- 
cernible offset or trend with respect to the luminosity 



the 



Karachentsev 



{2005) define 8 = max [log (M^/Df . )] + C to be 
amount a galaxy is acted on by its largest tidal disturber. As 
noted in that work, = corresponds to an object with a Keple- 
rian orbital period about MW/M31 equal to 1/H. For reference, 
Sgr dSph has = 5.6. 

5 If Fornax became a satellite of the MW relatively recently or, 
as is suggested by proper motion measurements, has a circular orbit 
about the MW, similarities in chemical evolution may be expected 
despite the current present day environmental differences. 



2003 



at all 



Woo et al 
when plot 



metallicity (LZR) relation of |Kirby et aL] ( | 2008 [). The off- 
set between dlrrs and d Sphs reported by jCrebel et al 

T 



2008) does not appear as strongly, it 
he updated metallicities for the 
galaxies we consider. Such systematic offsets are difficult 
to confirm and interpret in cases where biased metallic- 
ity indicators or small spatial coverage skew the average 
[Fc/H] values. WLM's agreement with the dSph LZR is 
consistent with studies looking at large samples of more 
distant galaxies from the SDSS survey - where only small 
deviations from the MZR are seen for galaxies of different 



envir onment or morphological class (Ellison et al. 2008 



20091. Environmental processing such as ram pressure 



stripping may not necessarily produce offsets in any one 
particular direction from the mass metallicity relation, 
however there is still likely a complex interplay between 
gas stripping, triggered star formation and dilution of 
the ISM during infall of a dwarf galaxy to a larger spi- 
ral. This may add to the observed scatt e r in the LZR 



(Skillman et al. 1996 
2W 



^ Boselli et al ||2008 
In addition varying levels of angu 



in the dwarfs may modulate SF efficiency ( Schroyen et al. 
20111, with the low angular momentum galaxies having 



Koleva et al. 



ar momentum 



more centrally concentrated gas which ma y be more effi- 
cient at cooling into the molecular phase ( Kuhlen et al. 
20111. 

The spread in metallicity, crrpe/H] i f° r t ne dwarf galax- 
ies studied here is nearly constant over 4 orders of magni- 
tude in luminosity, with a mean of (crpe/H]) — 0.38. The 
anticorrelation with luminosity found for the ultra faint 
dwarfs seems to saturate at this value for galaxies with 
luminosity L > 5 x 10 5 L Q in Figure [7J If the spread 
in [Fe/H] reflects the stochastic nature of enrichment 
events, then perhaps higher luminosity systems (where 
SF proceeds over several Gyr and many supernovae oc- 
cur) end up with a similar value of <7[Fe/H] which r eflects 



the mixing efficiency in the ISM (Argast et al. 2000 



Lower luminosity systems which have truncated SI 1 ' 
and therefore fewer enrichment events will unevenly dis- 
tribute metals through the galaxy in the short time they 
are forming stars leading to the higher dispersion [Fe/H]. 

5.2. Spatial Abundance Signatures of the Sample 

Past studies flWinnickl |2003[ |Koch et al.||2007 



lieuszik et al.|20 09 ; Kirby et al. 201l| have found bot 



Gul- 

nat 

and statistically significant negative radial [Fe/H] profiles 
for dSphs. However, in studies using linear CaT calibra- 
tions unbiased gradient estimates have proved difficult 
to compute. Even then small spatial coverage may mis s 
such a gradient. As pointed out by Kirby et al. (2011), 
their sample of stars in Fornax showed a significantly 
different gradient from the DART sample simply due to 
the small region of the galaxy sampled - which is rein- 
forced when looking at the inner regions of our Figure [9j 
Due to the strong RGB metallicity gradients shown by 
the dSphs, it is necessary to sample out to large radii in 
order to detect the low metallicity stars and build up a 
representative MDF. 

In this work we have selected literature data which 
show the least bias in terms of spatial coverage and 
[Fe/H] calibration to allow for a consistent differential 
comparison of the chemical trends with radius as traced 
by RGB stars. With the unbiased metallicity and spa- 
tial coverage in the sample we consider, there appears 
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to be a slight dichotomy in the abundance gradients of 
dlrrs and dSphs - with the gas rich rotating dlrrs show- 
ing significantly shallower radial gradients. While this 
could be due to differences in their total masses, fur- 
ther evidence supporting angular momentum as a driver 
of gradient strength comes when looking at the Local 
Group dEs: NGC 205, NGC 185, and NGC147. T he dis 
persiqn dom inated syste m NGC 205 (V/ a ~ 0.3; [Geha, 



et al. 



2006 ) was found by Koleva et al. ( 2009 ) to have a 
licity gradient comparable to the dSphs m our sam- 
> -0.14 dex r^ 1 ). By co ntrast NGC 



metal 

pler (d[Fe/H]/dr, 

147 a nd NGC 185 (V/a > 0.91, V/a > 0.65; |Geha et~al. 
2010| were found to have gradients as flat or flatter than 
the dlrrs we consi der here (d[Fe/H]/dr c ~ —0.02 dex r" 1 ; 
Geha et al.pOlOl ). 

'lb quantify this, in Figure [TTJ we plot the slope of the 
radial metallicity gradient versus V/a. V/a is not a sin- 
gular quantity, but changes with radius within a galaxy - 
therefore we show the range of V/a spanned by particular 
galaxies in this plot. The measured values for the dSphs 
(and the LMC) must be taken as coarse estimates, and 
should be interpreted with care. This is due to several 
factors: in the case of nearby dSphs with large angu- 
lar extents it can be difficult to ascertain whether the 
observed velocity gradients are intrinsic rotation or per- 
spective induced rotation, the uncertain inclination and 
axis of rotation for the dSphs may le ad to an underesti - 
mate of any true rotational velocity ( Lokas et al.||2010 ), 
and any tidal distortions to the dSphs make interpreta- 
tion of their dynamical state difficult. Nevertheless there 
is evidence for the lowest angular momentum systems to 
show steeper metallicity gradients in this literature sam- 
ple. 

It should be kept in mind that when comparing radial 
gradients (especially in systems with ongoing star for- 
mation) , ideally a tracer population of stars of the same 
age should be used. Indeed the change in abundance 
gradient as a function of population age offers a strong 



tool to stud y the chemical evolution (Cioni 2009} Vlajic 
et al.||2009|). For example the LMC shows a decrease 



in V/a with age, and the older stars are also more spa- 
tially extended. Therefore one would want to estimate 
the kinematics, chemistry and scale length of a single age 
stellar population in order to accurately test chemody- 
namic signatures - however this is extremely difficult to 
do with current data of Local Group dwarfs. The LMC 
may be a very complex case which should be analyzed 
with caution, however the similar metallicity gradients, 
dynamical state, mass and size of WLM, SMC, NGC 
147, and NGC 185 certainly offer a robust comparison 
sample to the dispersion dominated dSphs. While the 
dSphs may show on average older RGB stars, we note 
that WLM and Fornax have nearly identical age metal- 
licity relations (Fig. |Xo| ), yet show the most disparate 
radial abundance gradients. 

5.2.1. The Role of Environment 

The RGB gradients in Figure [9| and [Tl~1 for the gas rich, 
rotating dwarf galaxies are shallower than the low angu- 
lar momentum gas poor dSphs. However there are not as 
large differences between the metallicity gradient of the 
isolated WLM galaxy and the tidally disturbed Magel- 
lanic Clouds. If the metallicity gradient is modulated 
by environmental factors such as ram pressure or tidal 



stripping, this would be consistent with recent obser- 
vational and theoretical work suggesting that the Mag- 
ellanic Cloud s are on their first passage into the MW 



virial radius (Besla et al. 2007 Boylan-Kolchin et al 



20111) . Despite their close present day distance to the 
MW, their gas content and spatially extended SF may 
have been undisturbed, allowing chemical enrichment to 
proceed over the full body of the LMC and SMC for 10+ 
Gyr. In this environment driven scenario, the steeper 
metallicity gradients of dSphs would be due to their early 
infall times - with their gas content and chemical en- 
richment quenched as they were accreted by the MW. 
This could perhaps lead to a radially shrinking SF region 
(an d steeper metallicity gradient) within those dwarfs 
(i.e, |Mayer et al!|2007| |Mayer||2011 ) . 

In addition to ram pressure mediating metallicity gra- 
dients, it is not clear what factor tidal stripping of the 
stellar populations will play in pr e servin g or erasing a 



gradient. As shown by Sales et al. (2010) the kinematic 
and metallicity gradients in dSphs near the MW may be 
subject to tidally induced modifications depending on 
their initial strength and the orbital properties of the 
dSph. In most cases however, the metallicity gradient 
would remain, as the outer most metal poor stars would 
be unbound first. Therefore in situ formation of gradi- 
ents for dSphs could still be preserved (or made milder) in 
the prese nce of environmenta l processing from the MW 
(see also Koleva et al. 2011). The trends in Figure [Tl] 
may present a test tor merger or tid ally induced trans 
formations of dlrrs into dSphs (e.g, Mayer et al. 2001 



Klimentowski et al.|2009 Kazantzidis et al.|2011a|b 



as 

in such scenarios it could be difficult to produce obj ects 
with steep metallicity gradients from rotating progeni- 
tors which had flat radial metallicity profiles. Ram pres- 
sure or ad ditional baryonic effects may allo w for more 
flexibility ( |Mayer et al|2007|jK7avtsovj|2010[) on this is- 
sue (and on the evolution of the dwarfs' M/L ratios). 
Further modelling which jointly treats the chemical and 
kinematic evolution of dwarf galaxies in the MW poten- 
tial could comment directly. 

5.2.2. Internal Processes 



Simulations by Schroyen et al. (2011) showed that an- 
gular momentum plays a strong role in determining the 
radial metallicity profile of a dwarf galaxy. In their sim- 
ulations rotation produced a centrifugal barrier which in 
turn prevented gas from settling in the centre of galax- 
ies. This led to SF that occurred over the full extent of 
the dwarf at lower levels. This scenario naturally pro- 
duced smoother radial metallicity profiles, and extended 
star formation histories. These results were in qualita- 
tiv e agreement w i th the observations of dwarf ell ipticals 
by | Koleva et al. (2009), but a larger sample by Koleva 
etaT~(2011) from the literature has found weaker corre- 
latio ns with host galaxy pro perties. One simulation run 
from Schroyen et al. (2011) shows a metallicity gradi- 
ent of d[Fe/HJ/dr e ~ -0.12 dex rj 1 for the non-rotating 
dwarf, with the rotating dwarf (V/a ~ 1.8) having a gra- 
dient of d[Fe/H]/dr c 0.03 dex r" 1 . This is in excel- 
lent agreement with the average radial gradients we find 
forthedSphs (-0.13±0.01) and dlrrs (-0.03±0.01). To- 
gether with the dEs discussed above these observations 
would support a correlation between angular momentum 
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and radial metallicity gradient strength in Local Group 
dwarf galaxies. 

Another secular process for flatte ning metallicity gra- 
dients demonstrated in simulations flSellwood fc Binney 
20021 IRoskar et al.|[2008l |Stinson et al.p(M9j |Loebman 
et al.|201l[ ) is radial migration of stars. The migrations of 
the stellar populations are typically produced by global 
disk instabilities or transient spiral structures. These are 
thought to be common in MW sized galaxies, however 
it is unclear how ubiquitous such dynamica l instabili 



ties are in low mass, thickened dwarf galaxies ( Sotnikova 
fc Rodionov"! [20031 jSanchez-Janssen et al.||20l0| |Mayer 
201ip . Redistribution of locally enriched material within 
a galaxy is also possible due to the SF driven f o untain 
mechanism as shown by De Young & Heckman (1994), 
but aga in this may not be dom inant in low mass, thick 
dwarfs ( |Schroyen et al. 2011). If the dlrrs were high 
enough mass tor one or both of these processes to oper- 
ate, the dichotomy in metallicity gradients may then be 
due to total dynamical mass - however the similar masses 
of NGC 205, 185, and 147 would argue against this. 

In summary there could be several mechanisms at work 
to produce the observed radial [Fe/H] gradients. More 
data are needed to reveal whether environment, total 
mass, or angular momentum is the more fundamental 
parameter modulating metallicity gradients in dwarfs - 
which may be difficult as a correlation between V/a and 
environment may also be present. Efforts to incorporate 
metallicity profiles into dynamical simulations of DG evo- 
lution may also increase constraints on the contribution 
of environment. In all cases it is imperative to ascertain 
the dynamical and chemical profiles of the dwarfs out to 
large radii - as the value of V/a in particular will change 
dramatically from the inner to outer regions. 

6. SUMMARY 

In this paper we have presented [Fe/H] and age esti- 
mates for 126 RGB stars in the WLM dlrr galaxy. These 
complement the kinematic and structural study of WLM 
presented in our earlier papers, and represent some of the 
first spectroscopic abundances and velocities of individ- 
ual evolved stars in a truly isolated Local Group dwarf 
galaxy. For WLM and the literature data we have com- 
puted the metallicity distribution functions and charac- 
terized the bulk chemical abundance properties, as well 
as calculated the radial metallicity gradients, and age- 
metallicity relation for each galaxy as traced by the RGB 
stars. 

The key points from our study are: 

• The global metal abundances (mean, median) and 
MDF for the isolated dwarf galaxy WLM are sim- 
ilar to the SMC and the more luminous dSph For- 
nax - in expectation with the luminosity-metallicity 
relation for the Local Group. 

• WLM shows stars at as low metallicity as the dSph 
Fornax, however the dSphs on average show more 
extended metal poor tails, and are better fit by 
simple leaky box chemical evolution models than 
the dlrrs we consider. 

• The intrinsic spread in a galaxy's metallicity is con- 
stant over 4 orders of magnitude in luminosity, sug- 



gesting that the dispersion in [Fe/H] may saturate 
for dwarf galaxies with luminosities L> 5 x 10 5 L Q . 

• The dispersion in metallicity increases with age in 
WLM as is found in the LMC - with both galaxies 
also showing an increase in velocity dispersion in 
age (Paper II). 

• WLM, along with the SMC and LMC show radial 
[Fe/H] profiles that are statistically flatter than the 
dSphs. This along with the flat gradients for the 
rotating dEs NGC 147 and NGC 185 supports a di- 
chotomy in radial metallicty gradients which cor- 
relates with angular momentum for Local Group 
dwarfs. 

The strength of this study lies in the use of spatially 
extended samples of stars which have unbiased [Fe/H] 
indicators. Further measurements of stellar spectra at 
large radii may reveal as yet undiscovered velocity and 
metallicity gradients in dSphs or dlrrs which could help 
confirm the chemodynamic correlations suggested here. 
If borne out by further observations, such correlations 
between angular momentum and radial metallicity gra- 
dients could offer useful constraints for models of envi- 
ronmentally driven transformations of dlrrs into dSphs. 
Specific simulations tracking the angular momentum and 
radial chemical enrichment of an infalling dSph progeni- 
tor would be useful to improve the current understand- 
ing of how the classes of dwarf galaxies are connected. 
This coupled with additional spatially complete surveys 
of the chemistry in isolated dlrrs as well as dSphs will 
undoubtedly shed light on the magnitude of internal and 
environmental effects that shape the morphologies and 
chemistries of the Local Group dwarf populations. 
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Fig. 1. — Equivalent width as a function of magnitude relative to the horizontal branch for stars in the calibrating cluster s NGC 6791 
(red), Pal 14 (magenta), NGC 7078 (blue). Shown as the solid black lines is the [Fe/H] calibration of |Starkenburg et al.| |2010| l, illustrating 
the non-linearity clearly. Calibration lines proceed in constant [Fe/H] values from solar to —2.5 dex in U.5 dex steps, according to the given 
calibration. D ash ed red, magenta and blue lines are [Fe/H] values corresponding to the calibrating cluster's mean metallicity as taken from 
Harris (1996) and [Potter et al.| ( [2008| l (Pal 14). 
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,S2 dataset of Paper I are shown in blue, and DEIMOS spectra of the highest S/N quality in green. W 
tow good agreement, providing confidence in even the lowest S/N DEIMOS stars 
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Fig. 3. — Impact of [ Fe/Hl CaT calibr ation on the presence of spatially segregated subpopulations and gradients. Left panel shows 
the linear calibration of |Cole et a l. (2004) adopted in Paper I, compared to the non-linear calibration of Starkenburg ct al. (2010 1 in the 
right panel adopted for tms work. Probability shown in the lower panels are computed from a two dimensional, two-sided Kolmogorov- 
Smirnov test, and represents likelihood that the metal poor and metal rich populations are drawn from the same parent distribution in 
each calibration case. 
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FlG. 4. — Representation of systematic age errors on artificial input stars of various [Fe/H] and true age. Each oval shows an estimate of 
the combined systematic uncertainty in recovered age due to the effects of differential reddening, [o/Fe] variations, and AGB contamination 
for 10000 artificial star tests. Arrows show movement in the mean recovered age from the combined effects. 
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Fig. 5. — Panels show example of spatial and number biases in metallicity distributions of four dwarf galaxies in the sample. In addition 
to the full CDFs, the CDF within the inner 1.5r c is shown (dashed line), along with the mean and dispersion from 10000 iterations of 
random sampling of 31 stars from that inner region in each galaxy (blue solid and dashed). The small difference between the solid blue 
and black dashed lines indicates that the primary bias in comparing samples is consistent spatial coverage, especially when dealing with 
galaxies that show strong metallicity gradients. 
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Stars with r < 1 .5 r„ All Stars 




[Fe/H] [Fe/H] 

Fig. 6. — Differential MDFs for the Local Group dSphs and dlrrs considered in this paper, ordered by decreasing luminosity from top 
to bottom. Galaxy names are colour coded by the first author of the samples: red (Kirby), magenta (DART survey), orange (Magellanic 
Cloud surveys) (see Table 1). Left panel shows metallicity distributions in each galaxy where only stars within 1.5r c have been considered, 
right panels show full sample of stars for a given galaxy. Panels show normalized histogram only with respect to the particular sample (i.e., 
fraction of total stars for the right hand panels, or fraction of stars within 1.5r c for the left panels). Dotted blue, magenta, and red lines 
indicate the 10t h, 50th, and 90th percentile [Fe/H] values for each sample. Simple leaky box and pre-enriched chemical evolution models 
||Prantzosl2008) are overlaid on the full samples in green and cyan. 
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Fig. 7. — Intrinsic dispersion in [Fe/H] as a function of host galaxy luminosity for the sample of |Kirby et al.| ( [2010} (red circles) and other 
dwarf galaxies considered in this paper (filled black d ots). Solid lines connect values where the intrinsic dispersion has been measured on 
the same data in this work, and in |Kirby et al,| |2010) as a check on consistent methodology for removing the contribution of measurement 
uncertainty to the dispersion. Dotted lines connect the derived dispersion betwe en the different DA RT and Kirby et al. samples for Fornax, 
Sculptor and Sextans. Red dashed line shows the best fitting relation found by |Kirby et al,| ( |2010fr including the ultra faint dSphs. 
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Fig. 8. — Plot of [Fe/H] vs. geometrical radii in units of core radii for WLM and the Local Group dwarfs from literature. Galaxy names 
are colour coded by the first author of the samples: red (Kirby), magenta (DART survey), orange (Magellanic Cloud surveys) (see Table 
1). Dashed lines show weighted linear fits to the data. Solid lines show the running boxcar averages of Figure [9] overlaid. The rotating 
dlrrs appear to show statistically flatter abundance gradients than the dispersion dominated dSphs (weighted averages of —0.03 ± 0.01 vs. 
-0.13 ±0.01). 




Fig. 9. — Running boxcar averages of [Fe/H] as a function of radius in dlrrs (blue), and dSphs (red). Radial profiles have been smoothed 
by a factor of 2 for clarity, and normalized by the metallicity of the central regions of each dwarf galaxy. The blue and red shaded bands 
represent the associated 1 — a uncertainties on the running averages. 
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Fig. 10. — Age-metallicity relations for the LMC | |Cole et al"1|2004L SMC jCarrera et al.||2008a| l, Fornax dBattaglia et al.||2006| ), and 
WLM. Binned RGB field stars are shown as large black circles for a ll ga laxies, and in WLM the individual stars are shown as the small 
black dots. The oldest cluster in the LMC and WLM ( |Hodge et al.|1999| l are shown as the red triangles, and the A or B supergiant values 
as the blue squares. Overlaid in green and cyan are th e simp le leaky box and pre-enriched chemical evolution models from Figure [6] The 
HST SFH history solutions from the work of |Dolphin| ( |2000| l are overlaid as the magenta line for WLM. 
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Fig. 11. — The radial metallicity gradients plotted against the ratio of rotation to pressure support. The V/o lines in the x-dircction 
are not the errors, but show the range spanned by the galaxy over the radii where kinematic and metallicity information exist, and should 
be considered coarse estimates. In the dSphs it should be stressed that the observed rotation may not be intrinsic rotation, or may be 
overestimated due to uncertain kinematic axes or t idal feat ures, and should b e conside red coarse estimates o nly (see text). Da ta for the 
velocity info rmation comes from|Mateo et al.l |2008| (Leo I), [Koch et alj {2007} (Leo II), |Walker et a l. (2009b) (Fornax, Carina), iBattaglia l 
et al.H2008 t (Scu lptor) , |BattagIia et al.| ^lnBeitans), |GeTiaet^ 147, NGC 185), |OeTiaet~al.| ( |2006| l (NGC 205) , |Hams 

fc Zari tsky (2006} (SMC), and [van der Marel et al.| | |2002| (LMU). 
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TABLE 1 

Local Group Dwarf Galaxy Sample 



Galaxy 


Nstars 


Tmax 1 ft 


WLM 


180 


0.79 


LMC 


373, 59, 383 


0.78 


SMC 


349, 364 


1.14 


Fornax 


870 


1.20 


Sculptor 


629 


1.34 


Sextans 


180 


0.75 


Leo I 


825 


1.18 


Leo II 


256 


1.10 


Carina 


327 


1.06 



Reference 



Leaman et al. 



2 009]|2012[| 
; [PompOia 



Cole et al.|)|2UCl5|TfF6mpei'a et alT]||2008|;J 
<Jaxrera efa l. ( 2008ajj fFarisi et aL] j |2L)lb| r 
Battaglia et~al.| (|2U06| " 

Tolstoy ct al. (2004)'° 



Carrera et aL] ( |2008b[ ) 



Battaglia ct al. (201 
|Kirby et al.| I |2ul0| r 
|Kirby et"aT] |2UlUp 
Koch et al.fpJOlflp 



a Column shows what fraction of the tidal radius the outer most star in the spectroscopic sample extends 
to. 

b "DART sampl e" - original data from the se papers updated with additional observations and the [Fe/H] 
calibration from [Starkcnburg ct ah] ( |2010| ) . 



TABLE 2 
MDF Properties 



Galaxy 


10 th % 


50 th % 


90 th % 


V [Zq] 


[Fe/H] 


LMC 


-1.06 


-0.45 


-0.18 


0.430 
0.363 


— CO 

— 1.30 


SMC 


-1.53 


-1.05 


-0.64 


0.100 
0.085 


— oo 
-1.91 


WLM 


-1.74 


-1.24 


-0.75 


0.070 
0.064 


— oo 
-2.34 


Fornax 


-2.04 


-1.17 


-0.74 


0.093 
0.090 


— CO 

-2.66 


Leo I 


-1.84 


-1.42 


-1.11 


0.090 
0.037 


-co 
-2.18 


Sculptor 


-2.45 


-1.96 


-1.41 


0.014 
0.013 


-co 
-3.14 


Leo II 


-2.29 


-1.59 


-1.28 


0.036 
0.033 


-co 
-2.68 


Sextans 


-2.89 


-2.26 


-1.66 


0.007 
0.007 


-co 
< -5.0 


Carina 


-2.47 


-1.87 


-1.51 


0.019 
0.017 


— CO 

-3.02 



Note. — Effective yields in the first row for each galaxy 
represent the best fitting value from a Leaky box model, 
second row the effective yield in the pre-enriched model 
and initial [Fc/H]. 



